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Abstract 



We consider the problem of bounding away from the minimum 
value m taken by a polynomial P G Z [Xi , . . . , Xk] over the standard 
simplex A C M'"', assuming that m > 0. Recent algorithmic develop- 
ments in real algebraic geometry enable us to obtain a positive lower 
bound on m in terms of the dimension k, the degree d and the bitsize 
r of the coefficients of P. The bound is explicit, and obtained without 
any extra assumption on P, in contrast with previous results reported 
in the literature. 

1 Introduction 

1.1 Problem statement 

Let P G Z [Xi, . . . , Xk] be a multivariate polynomial of degree d taking only 
positive values on the /c— dimensional simplex 



Let r be an upper bound on the bitsize of the coefficients of P. Writing 



we consider the problem of finding an explicit bound mk^d,T depending only 
on k, d and r such that < mk^d,T < 



k 




i=l 



m = min P > 
A 



1.2 Previous work 



Several authors have worked on this subject. There are two main ap- 
proaches: Canny's gap theorem can be used, under non-degeneracy con- 
ditions ([C]); in |LSj . the authors use the Lojasiewicz inequahty, leading to 
a bound in the general case, but involving a universal constant. The method 
presented here gives an explicit bound, with no extra assumption on P. 



1.3 Univariate case 

We begin with the univariate case, which contains some basic ideas of the 
proof in the general case. This situation has already been studied in [BCR] . 
We present here a simpler proof, leading to a slightly better bound. 



Consider a univariate polynomial of degree d 



i=Q 

taking only positive values on the interval [0, 1]. Let r be a bound on the 
bitsize of its coefficients. 

The minimum m of P on [0, 1] occurs either at or 1, or at a point x* lying 
in the interior ]0, 1[. The first case is trivial, as P(0),P(1) € Z, so that m 
is clearly greater than 1. In the second case, P (x*) = 0, so that m is a 
root of the resultant R{Z) = Resr (P(T) - Z, P'{T)) E Z [Z]. The resultant 
R{Z) is the determinant of the matrix Syl(Z), where Syl(Z) is the following 
Sylvester matrix: 



/ 



ad 




{d-l)c 




ai 



ao — Z 



ad 



ai 



{d-l)ad-i 



\ 







ai uq — Z 






ai J 



d-1 



> d 



2 



d-l 



R{Z) = ^ rjZ* is thus a polynomial of degree d—1 in Z, whose coefficients 

i=0 

are controlled in the following fashion: 



Lemma 1.1. For all i £ {0, . . . ,d — 1}, we have 

d /d-l 



\ri\ < 



2^V(d+l) 



2^\/dTT- 1 



d-l-i 



Proof. Let (^i , . . . , ^d-i , Bi, . . . , B^i) denote the rows of the classical Sylvester 
matrix Syl(O). Then 



R{Z) = det (^1 + Zea+i, A^-i + Ze2d-i,Bi 



5 • • • ) ^aj ) 



B, 



where (ci, . . . , e2rf_i) is the canonical basis of M^"^^. Using the multilin- 

d~i 

earity of the determinant, we can write R{Z) = ^ r^Z*, where, for all 

i=0 

i £ {0, . . . ,d — 1}, ri is a sum of (^^7^) determinants of matrices built with: 

- i rows among the ej's 

- d — 1 — i rows among the Aj 's 

- the d rows Bi, . . . , B^. 

Hadamard's bound (see [BPRJ ) implies that, for all i: 



r,: < 



< 



d-l 



[{d + 1) (22- - 1)] 



d-l-i 



d{d + l){2d + l) 
6 



(22- - 1) 



d-l-i 



2^V(i+ 1 - 1 



d-l~i 



as claimed. 



d fd-l 

i 



2^VdTl - 1 



d-l-i 



□ 



Since the minimum m is a root of R{Z), Cauchy's bound finally implies 
the following theorem 

Theorem 1.2. Let P £ Z[T] be a univariate polynomial of degree d taking 
only positive values on the interval [0, 1] . Let r be an upper bound on the 
bitsize of the coefficients of P. Let m denote the minimum of P over [0, 1]. 
Then 

3d/2 



m > 



2(2d-l)r(^_^ 1^2d-l/2- 
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Proof. If m is attained at or 1, then the result is obvious. If not, m is a root 
of the resultant R{Z). Since R has at least one non-zero root {R{m) = 0), 
Cauchy's bound (see [BPRj ) implies 



i=0 
d-l 

i=0 



d /d-l 



2^V(i+ 1 - 1 



d-l-i 



1)3 



d r 



d-l 



from which the result follows easily. 



□ 



Remark 1.3. Our bound is slightly better than a recent one presented in 
\BCR^ . which was already almost sharp. Indeed, following \BCR^ . consider 
the polynomial = X'^ + (^2^X — l) . Here, t = 2k and the minimum 
of Pk satifies 

mk<Pk{2-'')=2-~''^l\ 
and thus decreases exponentially with d and r. 



2 Bound on the minimum of multivariate positive 
polynomial 

We now consider the multivariate case. 



2.1 Notation and problem statement 

The following notation will be useful: 

Notation 2.1. We write bit(n) for the bitsize of an integer n G N. 

Let P G Z [Xi, . . . , Xk] be a polynomial of degree d, and r a bound on 
the bitsize of its coefficients. Moreover, assume that 

m = minP > 0. 
A 

In order to find an explicit lower bound < mk,d,T < "i; we generalize 
the proof of the univariate case. We first show that, at the cost of slightly 
increasing the bitsize of the coefficients, we can assume that the minimum 
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is attained in the interior of the simplex. Obviously, there exists a face a of 
A, of dimension < s < fe, such that the minimum m is attained at a point 
of the interior of a (with its induced topology). In the following we consider 
such a face o", of minimal dimension s. 

Remark 2.2. If a is a vertex of A, then obviously m> 1. We now assume 
that s > 1. 



Denote by 



the vertices of A, and 



Vo = 

Vi = ei {l<i<k) 



Xo = l-EXi 
Xi=Xi {l<i<k) 

the associated barycentric coordinates. 

There exists a subset / = {io,. . . ,is} of {0, . . . , fc} such that the vertices 
of a are the vertices (l^)jgj. Let J = {0, . . . , k}\I. The face a is character- 
ized by: 

a = {x£ A |Vj G J, Xj{x) = 0}. 
We make the following substitutions in P: 

• li j e J and j > 0, replace the variable Xj by 

s 

• li j e J and j = 0, replace the variable Xi^ by 1 — ^ X^ 
We then obtain a polynomial ^TL [Xjj , . . . , Xj^] satisfying : 

min P = min P(j . 



Renaming the variables Xi^ into Y^, we obtain that P„ G Z [Yi, . . . , 1^] is a 
polynomial verifying: 

Lemma 2.3. The degree of P^ is hounded by d. 

Moreover, the bitsize of its coefficients is bounded by t^, where 

Ta = T + l + dhit{k). (2.1) 
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Proof. The degree of is clearly at most d. We now show the result 
concerning the bitsizes of the coefficients. 

The result is clear \i s = k. 

Assume that 1 < s < — 1. 
Since replacing Xj by does not change the bound on the bitsize of the coef- 

ficients, only the replacement of Yq by 1 — ^ has to be taken into account. 

i=l 



If 



then 



where 



and 



\o\<d 



\l\<d 



^ Ja(7i-/3i,.-,7fc-A)' 



= {/3 e N^+^ < d and e {1, . . . , s}, A < li}- 
Hence, we have 

m=p 

d 
p=0 

. {s + 



< 2' 



s 

d 



< 2^ X 2(s + 1) 



6 



and the conclusion follows. 



□ 



Since 



m = min P = min P. 



P(j G Z [li, . . . , Kj] achieves its minimum in the interior of o. Consequently, 
m is attained at a critical point of P^r-, i.e. a point x G such that the 
gradient of P^ is zero at x. We are thus interested in computing the values 
of P at the zeros of its gradient. We aim at giving a univariate reformulation 
of this problem, enabling us to use resultant methods. The following section 
introduces the necessary material. 

2.2 Rational univeiriate representation 

We first introduce the notion of Thom encoding: 

Definition 2.4. Let P G M[X] he a real univariate polynomial, x G M a 
real number and a G {0, 1, — a sign condition on the set Der(P) = 
(P, P', . . . , p(degJ')} of the derivatives of P. 

The sign condition a is a Thom encoding of x if a{P) = and 



We can now define a rational univariate representation as follows: 

Definition 2.5. An s— rational univariate representation u is an (s+S)— tuple 
of the form 



2. F and go are coprime, 

3. IT is a Thom encoding of a root t-,^ of F. 

Remark 2.6. IfteRisa root of F, then go{t) / 0. 

We now define the point associated to the rational univariate represen- 
tation: 




u=(P(T),5o(r),...,5.(r),7r) 



such that: 



1. F,go,...,gseR[T] 



7 



Definition 2.7. The point associated to u is defined by 

V5o(i7r)'"''5o(*7r)y ' 

Hence, a rational univariate representation gives rise to a point whose 
coordinates are rational fractions evaluated at a root of F. 

Let Q G M** be a nonnegative polynomial over M*, and 

Z{Q) = {xe R'\Q{x) = 0} 

be the set of real zeros of Q. We are interested in finding a point in each 
bounded connected component of Z{Q). This can be done by applying 
Algorithm 12.15 of [BPR] . which we recall here for convenience. 

Algorithm 2.8 (Bounded Algebraic Sampling). 

Require: A polynomial Q G ZlXi, . . . Xg], of degree bounded by dq, non- 
negative over M*. 

Ensure: A set U of rational univariate representations of the form 

{F{T),go{T),...,g,{T),7r), 

where the polynomials F, go, . . . , gg have integer coefficients, and such that 
the associated points meet every bounded connected component of Z{Q). 

We indicate the main ideas behind the algorithm, referring the reader to 
[BPR] for details: 

• replace Q by a deformation Def (Q, dq, () of degree bounded by dq + 2, 
where C is an infinitesimal, 

• consider the critical points Cr(Def (Q, dg, C)) of Dei{Q,dq,C) in the 
Xi-direction, 

• due to the properties of Def(Q, dq,C)i 

— Cr(Def (Q, dq, ()) has a finite number of points, 

— the quotient ring defined by the equations of Cr(Def((5, dq, ()) is 
a vector space of dimension at most 

{dq+2){dQ + l)''-\ 
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— its multiplication table can be easily computed, 



• find rational univariate representations of the points of Cr(Def (Q, dq, Q), 

• take their limits with respect to C, which define a finite set of points 
intersecting all the bounded connected components of Z{Q). 

The complexity analysis in [BPR] shows that, if dq is a bound on the degree 
of Q and tq a bound on the bitsize of its coefficients, then: 

1. The degrees of the polynomials F, qq, . . . , are bounded by 

{dq + 2){dq + l)'-' 

2. The bitsize of their coefficients is bounded by 

{dq + 2){dq + 1)''-Hkdq + 2) (r' + 2bit(A;dQ + 3) + 3/i + bit(A;)) , 
where 

r' = sup [tq, bit(2A;)] + 2 bit [k{dq + 2)] + 1 



H = bit 



{dq + 2){dq + l 



2.3 The bound 

Recall that G Z [Yi, . . . , Y^] achieves its minimum in the interior of a. 
Consequently, this minimum is attained at a critical point of P^, i.e. a point 
X G at which the gradient of P^ is zero. Consider the set of critical points 

^ = |.,r|||(.) = ... = §|w = o}. 

Note that Pa is constant on each connected component of Z. So, we aim at 
computing a set U of rational univariate representations u whose associated 
points Xu meet every connected component of Z, together with the values 

Remark 2.9. When Z has a finite number of points, then Grohner basis 
techniques can be used to obtain rational univariate representations of these 
points (see fB^). The method we present hereafter, based on Alaorithm \2.8[ 
computes a point in every connected component of Z even if Z is infinite. 
Moreover Alaorithm \2.8\ makes it possible to control the degree and the bitsize 
of the coefficients of the output, in contrast with Grobner basis methods. 
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It is easy to see that if C is a connected component of Z containing 

o 

a minimizer of Pa in o", then C C o" by minimahty of the dimension s of 
a. In particular, C is bounded. Algorithm 12.81 then gives a set of rational 
univariate representations of the form 

^x = (F(T),9o(r),<7i(T),...,<7s(T),7r), 

whose associated points meet every bounded connected component of Z. In 
particular, they meet every connected component of Z containing a mini- 
mizer of P(j in fj. 

Lemma 2.10. The degree of the polynomials F, go, . . . , gg is bounded by d^, 
where 

du = 2d{2d-lf-\ 
Moreover, the bitsize of their coefficients is bounded by 

Tu = du{2kd-2k + 2) [r' + 2 bit(2/c(i - 2/;: + 3) + 3 bit(d„) + bit(A:)] , 

where 

t' = 2T + (2d + 2)bit(/c) + (A: + 3)bit((i) + 5. 
Proof. Let Q denote the polynomial 



i=l 



Clearly, its degree is bounded by dq = 2d — 2. Moreover, we can bound the 

bitsize of its coefficients as follows. 

If 

P.= Y. 

\l\<d 



then 



b\<d 

where 

a<-f 
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Write Q = J2 dsY^ . Since Q = J2 { "^^rr i its coefficients are bounded 



as follows: 



5eN' i=l 
\S\<d-2 



dYi 



am" 
a<'Y 

Hence, the bitsize of the coefficients of Q is bounded by tq, where 

TQ = 2r„ + (/c + 2) h\t{d) + bit(A;) = 2t + {2d + 1) bit(A;) + {k + 2) h\t{d) + 2, 



where the last equality follows from equation (12. ip . The result now follows 
from the complexity analysis of Algorithm 12.81 □ 

Let Pu{T) = goiTfP. f • • • , ^) • We have: 

\9o{T) 9o{T)J 

Lemma 2.11. The degree of Pu is bounded by dp^u = dud. 
The bitsize of its coefficients is bounded by Tp^u, where 

rp,u = d[Tu + bit((i„ + 1)] + r + dbit(A:) + d + k + l. 

Proof. The result about the degree is clear from the previous lemma. 

The bound on the bitsize of the coefficients is obtained by substitution, using 

Proposition 8.11 of [BPR] . □ 

The minimum m of P^ over a is attained at a point x € a contained 
in a connected component of Z included in the ball -6(0, 1). Since P^ is 
constant on such a component, m is also attained at some point Xu as- 
sociated to an already computed rational univariate representation u = 
(F(T),5o(T),5i(T),...,ff,(r),7r). 

Since tj^ is a root of F, the minimum m = P^ (xu) is a root of the 
resultant 

R{Z) = ReST (Pu{T) - go{TfZ,F{T)) . 



Example 2.12. We consider here the following easy example (Berg polyno- 
mial, see Example 37 in fSB/): 



B := x^y'^ix^ + y'^ - 1) + 1. 
It is easy to show that B is positive on A. We now compute its minimum. 
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On the three vertices of A, we have B = 1. 

On the faces {x = 0} and {y = 0}, we have B = 1. 

Consider the face {x + y = 1}. Replacing x by I — y leads to consider 
the (univariate) polynomial 

Since B'^^^^^-^^^ = 6y^(y — l)^(2y — 1), the minimum of B^^_^y^ij is 
31/32, attained at y = 1/2. 

We now compute the values of B at its critical points contained in the 
interior of A. It is easy to show that those points {x,y) satisfy 

2x^ + = 1 
+ 2y2 = 1. 

Note that this easily implies that there is only a finite number of such 
critical points. A rational univariate representation of this set can then 
be computed (using for example Salsa software, see fSa]/): 

F = (3r2 - l)(r2 -3) 

50 = T(3T2 - 5) 

51 = T2 + 1 

52 = 2(t2 - 1). 

The resultant R{Z) is equal to 

mZ) = R«. (*(T)«B (ff}. ff}) - Z,„(Tf, F(T)) 
= 2^^3^(27Z - 2Qf. 

The only root 26/27 < min(l, 31/32) is thus the minimum of B over 
A, corresponding to the root \/3 of F , and giving the minimizer 

gi(V3) 52 (^) \ ^ M ^\ 
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In order to obtain a lower bound on the minimum depending only on k, d 
and T, one needs to bound the roots of R{Z). This can be done by controlling 
the size of the coefficients of R{Z) and then using Cauchy's bound. Write 



F{T) 



du 



1=0 



dp,u 

PuiT) = ^aiT' 



1=0 



and 



dud 



dp,t 



90 



i=0 



i=0 



Lemma 2.13. The polynomial goiT^ has degree hounded by d^d = dp^u, 
and the bitsize of its coefficients is bounded by d{Tu + bit((i„ + 1)) < Tp^u- 

Proof. The degree of go{TY is clearly less than d^d = dp^^. We now show 
the bound on the bitsize of its coefficients. Recall that the degree of g^ is 
bounded by d^and that the bitsize of its coefficients is less thanxu. When 
multiplying a univariate polynomial / by go, the increase in the bitsize of 
the coefficients is at most r„ + bit((i„ + 1). Indeed, the coefficients of fg^ 
are sums of at most + 1) products of a coefficient of / by a coefficient of 
5o- The conclusion follows easily. □ 

The resultant R{Z) is the determinant of the matrix Syl(Z), where 
Syl(Z) is the following Sylvester matrix: 



(ad 



p,u ~ ^dp,u^ 




fdu 





ao — hoZ 



adp^^-hdp^^Z ■■■ 
fo 



fdu 



\ 





ao - hoZ 






fo J 
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du 

R{Z) = ^ rjZ* is a polynomial of degree < du in Z, whose coefficients are 
controlled in the following fashion: 



Lemma 2.14. For all i £ {0, ... , du}, 

r -\du 



ri < 



dp.u 



Proof. We proceed as in the proof of lemma ll.li 
Let (with obvious notation) 

{Ai + ZBi, . . .,Ad^ + ZBd^,Ci, . . . ,Cdp J 

be the rows of the Sylvester matrix Syl(Z). Using the multilinearity of the 
determinant, we can write R{Z) = ^ r^Z*, where, for all i £ {0, . . . , du}, r^ 

i=l 

is a sum of C^") determinants of matrices built with: 

- i rows among the -Bj's 

- du — i rows among the Aj^s 

- the dp^u rows Ci, . . . , Crf^^. 

Hadamard's bound (see |BPRj ) implies that, for all i: 
^du 



\ri\ < 



< 



^ y [(dp,„ + 1) (22--- - i)]'V[(^- + 1) (2'^" - 

\ r / 1 '^ti r / 



as claimed. 



□ 



Since the minimum m is a root of R{Z), Cauchy's bound finally gives 
the estimate we were looking for. 

Let U he a set of rational univariate representations whose associated 
points meet every bounded connected component of 



X G 



dPa. 



dPg 



ix) = 



for each face cr of A. 
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Also, let du (resp. r^) be a bound on the degree (resp. the bitsize 
of the coefficients) of the polynomials occuring in the rational univariate 
representations of U, and dp^u (resp. rp^u) be a bound on the degree (resp. 
the bitsize of the coefficients) of the polynomial Pu- Then: 

Theorem 2.15. 

1 

Proof. Since R has at least one non-zero root {R{m) =0), we can write 

p 

R{z) = Y,nZ\ 

i=q 

with q < p < du and VqVp ^ 0. 



Cauchy's bound then implies: 

du 



1 

- < El 

i=0 
du 

< E 

i=0 



du r 



2"^-V^P,« + l 2^■-^/du + l 



dp. 



< 2'^- \rp'- ^dp^u + 1] '^^ [2^- ^du + 1 



< 



du r 



as announced. 

Using Lemmas 12.101 and 12.111 Theorem 12.151 immediately implies: 



□ 



Theorem 2.16. Let P £ 1^ [Xi, . . . ,Xk] be a polynomial of degree d, t a 
bound on the bitsize of its coefficients and m = rmn P the minimum of P 

over the simplex A. Assume that m > 0. 
Let 

D = 2d{2d - l)''-^ , 

p = D{2kd -2k + 2) [r' + 2 hit{2kd - 2A; + 3) + 3 hit{D) + hit{k)] , 
p' = d[p + bit(Z) + 1)] + r + dbit(A:) +d + k + l, 
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where 

t' = 2T + (2(i + 2)bit(A:) + (A: + 3)bit((i) + 5. 

Then: 

m > mk4,r, 

where 

_ 1 

" [2p'+iVZdTT]''[2p/D+T]'^' 

Remark 2.17. We now give a more compact bound, derived from the last 
corollary. This will enable us to compare our results to those of de Loera 
and Santos (ILSj/ ) and Canny (fDil)- 
The following estimates 

D + l< 2^d^ 

2k+fd'^+2k - 2^ + + ^) + + ^) ^^*(^) + 6fc + 9 

lead to the bound: 

1 < (^2'^ '^2''"^^ <i''+^k 22''+^ dl'+^k^ 1^2''+^ d''+^k^2''+5d''+'-k^ ^2 2) 

/ri l-LiSy , t/ie authors give the following estimate: 

^ < (2-)^^'^"^' 2^^"="^ (t) 



rrikd- 



where B denotes a bound on max((i+l, fc+l) andc is an (unknown) universal 
constant. Note that bound ()2.2p implies that 

1 ^ /2rN2''-+»d'=+iA:22''-+''d'=+2fc2 



giving an explicit version of estimate ([f]) . 

Moreover, a direct application of Canny 's theorem (f^). under a non- 
degeneracy assumption on the following polynomial system (with unknowns 

' P^{Yi, . . . ,Ys) = m 

^(Y, Y)-...-^(Y Y)-0 
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leads to the estimate 



In the general case, (|2.2p implies 



1 , J J s 9'=+5i.2 jfe+l 



rrikA. 







where the main difference with ^ is the presence of the exponent 2^~^^ . This 
essentially comes from doubling the degree of the system ([S|) by replacing the 
equations 

dPg _ _ dP^ 
dYi "' dYs 

by the following single one 

2 

1=1 

in order to cover the degenerate cases as well. 
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